High- Temperature Atomic Superfluidity in Lattice Boson-Fermion Mixtures 
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We consider atomic Bose-Fermi mixtures in optical lattices and study the superfluidity of fermionic atoms due 
to s-wave pairing induced by boson-fermion interactions. We prove that the induced fermion-fermion coupling 
is always attractive if the boson-boson on site interaction is repulsive, and predict the existence of an enhanced 
BEC-BCS crossover as the strength of the lattice potential is varied. We show that for direct on-site fermion- 
fermion repulsion, the induced attraction can give rise to superfluidity via s-wave pairing, at striking variance 
with the case of pure systems of fermionic atoms with direct repulsive interactions. 
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The recent spectacular acceleration in the experimental ma- 
nipulation I; jj, y, y, |4J] and theoretical understanding |5|, |6[ of 
systems of neutral atoms in optical lattices is leading to new 
and far reaching possibilities in the study of complex systems 
of condensed matter physics. Besides opening the way to the 
controlled simulation and experimental testing of models of 
strongly correlated systems OLD* sucn as > e -g-' high-T c su- 
perconductors and Hall systems, atomic physics in optical lat- 
tices hints at the possibility of discovering and probing new 
quantum phases of matter. Among others, very interesting re- 
cent studies in this direction are concerned with lattice boson- 
fermion mixtures HUES EH El- 

Optical lattices are stable periodic arrays of microscopic 
potentials created by the interference patterns of intersecting 
laser beams. Atoms can be confined to different lattice sites, 
and by varying the amplitude of the periodic potential, it is 
possible to tune the interatomic interactions at will. There- 
fore, optical lattices provide an ideal tool to reach strong cou- 
pling regimes even in the dilute limit Q] ■ The theory of neu- 
tral bosonic atoms in optical lattices has been developed 
by assuming that the atoms are confined to the lowest Wan- 
nier band of the periodic potential. It can then be shown 
d that the system is effectively described by a single-band 
Bose-Hubbard model Hamiltonian 11311 . In such a model the 
superfluid-insulator transition is predicted to occur when the 
on-site boson-boson interaction energy becomes comparable 
to some multiple of the hopping energy between adjacent lat- 
tice sites. This situation can be experimentally achieved by 
increasing the strength of the lattice potential, which results in 
a strong suppression of the kinetic (hopping) energy term. In 
this way, the superfluid-Mott-insulator quantum phase transi- 
tion has been realized in a series of beautiful experiments by 
loading an ultracold atomic Bose-Einstein condensate in an 
optical lattice d. The theory of dilute mixtures of interact- 
ing bosonic and fermionic neutral atoms subject to an optical 
lattice has been recently developed by deriving an effective 
single-band Bose-Fermi Hubbard (BFH) Hamiltonian from 
the underlying microscopic many-body dynamics (§]. The 
zero-temperature phase diagram of the BFH model has been 
analyzed both in the homogeneous and inhomogeneous cases 

SUES. EH El 

These recent studies have shown that Bose- 
Fermi mixtures are very fundamental systems of condensed 
matter exhibiting complex phase diagrams and a potentially 



very rich physics that is just beginning to be unravelled. 



The quest for superfluidity of fermionic atoms in dilute, de- 
generate quantum gases is the current o bject of widespread 
theoretical and experimental efforts 0, 11 El- In the 
present work we propose a new mechanism for fermion pair- 
ing in atomic systems. We consider mixtures of bosonic and 
fermionic atoms in optical lattices. The fermions are assumed 
to belong to the same atomic species, but trapped in different 
spin states, so that they can differ by angular momentum or 
generalized spin a (a = {f , J,}). This situation can be easily 
achieved, for instance, by all-optical trapping. The mixture 
is then characterized by independent bosonic and fermionic 
hoppings and by direct on-site boson-boson, boson-fermion, 
and fermion-fermion couplings. For such a system we deter- 
mine the structure of the induced fermion-fermion coupling 
due to the presence of the bosons, in analogy with recent 
and important studies on the effects of induced interactions 
in boson-fermion mixtures in free space 11611 . We show that 
the regime of static linear response is an excellent approxi- 
mation when the bosons are fully superfluid (weak optical lat- 
tice potential for the bosons) and, provided that the bosons are 
endowed with repulsive self-interactions, we find that the in- 
duced on-site fermion-fermion interaction is always attractive, 
thus leading to local pairing, irrespective of the sign of the 
boson-fermion interaction. In the case of a direct, repulsive 
on-site fermion-fermion interaction we show that the addition 
of the induced, attractive fermion-fermion interaction leads to 
an enhanced crossover regime between extended BCS pairing 
and Bose-Einstein condensation of local molecules, and al- 
lows for the occurrence of s-wave superfluidity of fermionic 
atoms. This is at striking variance with the case of pure 
Fermi systems with repulsive interactions, for which s-wave 
pairing is impossible and superfluidity might be achieved, at 
much lower temperatures, via higher-order mechanisms such 
as anisotropic <i-wave pairing. 



The single-band BFH system in an optical lattice with 
nearest-neighbor hopping and on-site direct boson-boson, 
boson-fermion, and fermion-fermion interaction is described 
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by the Hamiltonian | 



H 
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where Jb and Jf are, respectively the boson and fermion 
hopping amplitudes between adjacent lattice sites, and 
Ubb, Ubf, and Upf are, respectively the direct boson- 
boson, boson-fermion, and fermion-fermion on-site interac- 
tion strengths. The operators &i are the on-site bosonic anni- 
hilation operators, and Ci >a are the on-site fermionic annihila- 
tion operators for states with spin a. The bosonic on-site oc- 
cupation number is nf = a\di, while hi^ = c' icr Ci il7 denotes 
the fermionic on-site occupation number for states with spin 
a, and, finally, nf = fti^ + n^i denotes the total fermionic 
on-site occupation number. The bosonic and fermionic lattice 
potentials Vb and Vf are produced by a standing laser wave: 
Vb,f(x, y, z) = Vb,f (sin 2 kx + sin 2 ky + sin 2 kz), where 
k = tt/o, is the wave vector of the laser field, a denotes the 
lattice spacing, and the lattice stregths V b ,f that are felt by 
the bosons and the fermions depend on the intensity of the 
laser light and the detuning of the laser fequency against the 
respective atomic frequencies. Scaling the tunable lattice am- 
plitudes Vb,f with the bosonic and fermionic recoil energies, 
respectively E§ = (nh) 2 /2a 2 m B and E F = (irh) 2 /2a 2 m,F, 
one can define the dimensionless bosonic and fermionic lat- 
tice amplitudes n B ,F = Vb,f / E§ f- Gi yen the wavelength 
Al of the laser field, and the internal-state wavelengths As, 
Xf for bosons and fermions, one has np/r/B = ((Af/Ab) 4 - 
(T F A\ B E§))/(r B AX F E§), where AA B , F = Xb.f - A L 
are the boson and fermion detunings from the laser field, 
and the ratio of the natural line widths Tp/Vs is usually 
of order one. By properly adjusting the detunings (blue or 
red) with the gauging of the laser field wavelength Al, one 
can realize lattices whose depths can either be very different 
or nearly equal for the bosons and the fermions. The hop- 
ping and on site energy amplitudes exhibit a characteristic 
behavior as a function of the parameters of the optical lat- 
tice: J b ,f = E§ F (2/ y fjr)ri B /4 F exp(-2'q 1 J 2 F ); U B b,ff = 

E B F a,BB ,FFk^J% / irrfyp , where cibb and off are the 
direct boson-boson and fermion-fermion s-wave scattering 

lengths, and Ubf = EfaBFk(A/ ^){l+m F /m B ){n~ F 1 ^ 2 + 

rjg 1 ^ 2 )^ 3 ^ 2 , where asF is the direct boson-fermion s-wave 
scattering length. The ratio of the hopping amplitudes reads 

Jb/Jf = [(E§vT)/(E^vT)} ■ ex P [-2^- ^f¥)]- 

The bosonic and fermionic dispersion relations in a 3-D op- 
tical lattice with nearest-neighbor hopping take the following 
form: 



£ B,iKq) = 4J b .f 
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hopping amplitudes. The bosonic energy spectrum in the Bo- 
goliubov approximation reads 



^s(q) = V e B{q) [es(q) + 2n B U BB ] ■ 



(3) 



In the low momentum limit for the bosons, we recover 
the standard quadratic single particle spectrum es(q) as 
3 JBa 2 q 2 , where q 2 = \q\ 2 = q 2 . + q 2 + q 2 . Correspondingly, 
the phononic Bogoliubov excitation spectrum in the long 
wavelength regime reads ws(q) ~ \/ '3 J ' B{o,q) 2 2nBV bb = 
esq, where cb — ^/GubJeUbbo. is the velocity of sound. 

Let us next consider the momentum and frequency depen- 
dent bosonic dynamic response function xb(u, q): 



Xn(o>,q) = 



2n B e B {q) 
^ 2 - w|(<?) 



(4) 



Regarding the effects induced on the fermions, it is possi- 
ble to show that the bosonic dynamic response function can 
be replaced by its static limit approximation. Let us focus 
on the fermion-fermion interaction Upp(u), q) induced by the 
bosons. If the frequency contribution is not neglected, the gap 
function A(w,q) for the fermions depends on frequency av- 
erages over the off-diagonal fermionic propagator Gi2(u>, q). 



A(w,q) = J 'u$#(w' -w.q' - q)G 12 (w',q')dw<V . (5) 

The maximum scale of integration for the fermions is fixed 
by the Fermi energy: hio max as Ep. The condition w 2 riaa . <g; 
u B (q)\ q= q F can be satisfied for widely different scenarios of 
fermion filling factors. We will consider half-band fermion 
filling np as 1. This amounts to qpa as ir/2 and Ep as 6Jp. 
We then have T?uj 2 max as 36J F , and uj 2 B (q)\ q=qF as 36J| + 
12 JbkbUbb, so that the static limit condition reads 



3J F < 3J| + JbtibUbb 



(6) 



where Jb and Jp are, respectively, the bosonic and fermionic 



This condition can be satisfied in several ways, depending on 
the type of boson-fermion mixture that one actually consid- 
ers, and on the lattice strenghts and hopping amplitude ra- 
tios that can be realized in experiments. For boson-fermion 
mixtures with almost equal atomic masses, like 6 Li- 7 Li, and 
large red detuning of the laser field from the typical D lines 
of the Lithium isotopes, one has Jp ~ Jb, which implies a 
static limit condition of the form 6(Jf — Jb) "C nsUsB- 
This condition is well satisfied at the onset of bosonic super- 
fluidity Ubb — 10 Jb and bosonic filling factors n B > L 
Alternatively, the static limit condition can be satisfied by 
realizing small detunings such that rjp > 2t]b', this yields 
in general Jb ^ Jf, which again allows inequality l|6} to 
hold. In this second instance, the bosonic lattice is usu- 
ally some times weaker than the fermion lattice, and the 
bosons are in the deep superfluid regime. This setting is 
suited to treat mixtures irrespective of the atomic mass dif- 
ferences, and is thus the one we will consider in the follow- 
ing when studying the degenerate mixture of 40 K and 87 Rb. 
We will study quantitatively the case of half-band fermion 
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filling ray w 1, because this is a situation that can be eas- 
ily realized experimentally, it provides the optimal scenario 
for the highest attainable critical temperatures, and finally al- 
lows a direct comparison with the work of Hofstetter et al. 
on the high-temperature superfluidity of pure fermionic sys- 
tems in optical lattices 0]. Taking the static bosonic re- 
sponse function Xb(0, q) = —2nBes(9)/w^((7), the result- 
ing induced fermion-fermion interaction (due to the bosonic 
density-density fluctuations caused by the boson-fermion cou- 
pling U B f) reads Up$(0,q) = Xb(0><1)U bf . To discuss 
s-wave fermion pairing and evaluate the critical temperature, 
one must compute the average induced interaction Upp over 
the Fermi surface Jl6ll . Provided that the bosonic lattice is not 
extremely weak (i.e. r\B not -c 1), contributions from nearest- 
neighbor induced interactions are negligible in first approx- 
imation, and one finds U™ F = — U bf /Ubb- We see that, 
irrespective of the sign of the boson-fermion s-wave scatter- 
ing length cibf, and thus of the boson-fermion on-site cou- 
pling U BF-, the induced fermion-fermion interaction is always 
attractive as long as the bosons are endowed with a positive 
s-wave scattering length a^s > and thus with a repulsive 
on-site interaction Ubb > 0. This is the case we will always 
consider in the following. For distinguishable fermions, such 
as ensembles of spin unpolarized neutral atoms, the total on- 
site fermion-fermion coupling Up° F is then 

u 2 

jrtot rrdir | jjind jjdir BF /n\ 

u FF = u FF ~r u FF — u FF ~ 77 ■ \' ) 

Ubb 

Thus, regarding interatomic pairing, fermions in a lattice 
Bose-Fermi mixture behave as interacting particles in a pure 
lattice fermionic system, but for the crucial difference of ac- 
quiring a dressed attraction that modifies the direct fermion- 
fermion interaction through the boson-fermion coupling. This 
has relevant consequences on the possible pairing mecha- 
nisms in the system, in view of the relative sign between 
the direct and the induced part of the coupling strength. Be- 
cause the sign of U FF is determined solely by the sign of the 
fermion-fermion s-wave scattering length <Xff> we need to 
treat separately the two possible instances app < of direct 
on-site attraction, and app > of direct on-site repulsion. 

I) Fermions with direct on-site attraction - In this case 
the fermions can undergo a transition to s-wave superfluid- 
ity at a higher critical temperature T c than the analogous pure 
Fermi case, as the direct and induced interactions are both 
attractive and add together (See Eq. (Q). All the consider- 
ations valid for pure fermionic systems fTvtl hold as well 
in the case of a Bose-Fermi mixture, with the crucial dif- 
ference that the thermodynamic properties will now depend 
on the total fermion-fermion interaction Up° F rather than the 
direct one U FF . Starting with a lattice of low or interme- 
diate depth, the BCS picture for s-wave pairing holds, pre- 
dicting a critical temperature for the transition with a scal- 
ing behavior ksT c ~ 6Jp exp (— 7Jp/\Up°p\). As the depth 
of the lattice is increased, the on-site interaction becomes at 
first comparable and then finally dominates over the tunnel- 
ing amplitude; then the fermionic atoms form localized on- 
site bosonic molecules that can undergo a Bose-Einstein con- 
densation (BEC) into the superfluid state at a transition tem- 
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FIG. 1: Critical temperature for the SF transition of 6 Li atoms as 
a function of the fermion lattice strength (dotted line), and of 6 Li 
mixed with bosonic 7 Li (continuous line). 



perature ksT c ~ 6Jp/\Up° F \. In this situation the reduced 
ability of the pairs to move around in the lattice leads to a 
net decrease of the critical temperature. The maximum of the 
critical temperature T™ ax is then achieved at the region of 
crossover between the BCS and BEC regimes, |J7j?J.| — 10 Jf 
I17lll8ll . In Fig.^we show the behavior of the critical temper- 
ature for the superfluid (SF) transition of 6 Li atoms as a func- 
tion of the strength of the optical lattice, at np ~ 1, When 
trapped in the state | |> (4))> the fermionic Lithium atoms 
are endowed with a very large and negative scattering length 
clff — — 2500ao- This is the most favorable case of direct 
fermion-fermion attraction known so far, and considered by 
Hofstetter et al. 0]. Their result is compared in Fig. \l\ with 
the case of a 6 Li- 7 Li mixture, for which asF — 38ao QjJ. 
We see that the presence of the bosons further enhances the 
maximum critical temperature attainable for the SF transition 
of the fermionic atoms. For a lattice spacing a = 10 4 ao the 
maximum critical temperature T" Lax = 0.11/iA' for pure 6 Li; 
for the mixture 6 Li- 7 Li we have T™ ax = O.UfiK. 

II) Fermions with direct on-site repulsion - In the case of di- 
rect on-site repulsive fermion-fermion interactions, U FF > 0, 
the picture departs more radically from the analogous pure 
Fermi system. The direct repulsion and the induced attrac- 
tion that determine the total interaction strength are now in 
competition, and the net total sign depends on their relative 
importance. In fact, because U FF and Ubb are both pos- 
itive, we can evaluate the total fermion-fermion interaction 
as Up°p = (UbbUpp — U bf )/Ubb- In several important 
instances of mixtures of alkali atoms of current experimen- 
tal interest, with direct fermion-fermion repulsion, the mea- 
sured value of the boson-fermion s-wave scattering length 
clbf turns out to be much larger than both the boson-boson 
and fermion-fermion s-wave scattering lengths a BB and ciff, 
with the largest difference attained in degenerate mixtures of 
fermionic 40 K and bosonic 87 Rb. Namely, in this case one 
has that a FF = 104.8 ± 0.4ao (ao being the Bohr radius) in 
the singlet state, and a FF = 174 ± 7ao in the triplet state 
j2(|, while clbb = 100.2 ± 1.8a , and ag F — -410 ± 80a 
12111 . A similar situation occurs for a mixture of fermionic 
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FIG. 2: Critical temperature for the SF transition of 40 K atoms mixed 
with 87 Rb as a function of the fermion lattice strength without a Fes- 
hbach resonance (dotted line) and with a Feshbach resonance at 725 
G (continuous line). 

6 Li with direct repulsion and bosonic 7 Li, where one has 
cibb — apF — 5ao, and clef — 38ao. In these important 
instances the direct fermion-fermion repulsion is thus dom- 
inated by the induced fermion-fermion attraction, and high- 
temperature superfluidity of the fermionic atoms occurs by 
s-wave pairing. This is at striking variance with the case of 
systems of pure fermionic atoms with direct on-site repulsion: 
for such systems s-wave pairing is obviously always forbid- 
den, and superfluidity can be achieved via different mecha- 
nisms of higher order, at lower critical temperatures, such as 
anisotropic d-wave pairing. Therefore, mechanism II) real- 
izes the first atomic analogue of electron paring in crystal lat- 
tices. A further dramatic increase in the transition tempera- 
ture T c can be obtained by combining the boson-induced in- 
teraction with the experimental manipulation of the s-wave 
boson-fermion scattering length via Feshbach resonances. In 



fact, the existence has been recently predicted l22ll of a Fes- 
hbach resonance at an applied magnetic field B ~ 725 G in 
the degenerate mixture of 40 K and 87 Rb, leading to an en- 
hancement of the boson-fermion s-wave scattering length to 
the very high value gbf — — 687ao- In Fig.[2]we show the 
behavior of the critical temperature for the superfluid transi- 
tion of 40 K atoms induced by the presence of the 87 Rb atoms. 
Assuming the isotopes trapped at the Z?-line values of current 
experiments (795nm for rubidium, 767nm for potassium), and 
a blue-detuned laser at 764nm, we have rj F ~ Aj]b, so that in 
the range of values of r\p ~ 7 around the maximum criti- 
cal temperature, we have t\b — 2 (fully superfluid bosons) 
and J| ~ 20 J|. (static regime holds). We compare results 
for the case with (continuous curve) and without the presence 
of a Feshbach resonance (dotted curve). For a lattice spac- 
ing a = 10 4 ao the first case yields T" lax = 0.05/j.K, while 
without resonance we have T™ ax — Q.OifiK. In conclusion, 
we have shown that for mixtures of bosonic and fermionic 
neutral atoms in optical lattices, the fermionic atoms acquire 
an induced on-site attraction mediated by the boson-fermion 
interaction, irrespective of the sign of the latter. This fact al- 
lows the fermions to undergo a high-temperature superfluid 
transition via total s-wave pairing in both instances of direct 
attractive and repulsive on-site fermion-fermion interactions. 
Recent progresses in the manipulation of mixtures of bosonic 
and fermionic atoms in optical lattices 0] indicate that these 
predictions may be subject to experimental verification in the 
near future. We warmly thank Dr. K. Bongs for providing us 
with up-to-date values on forthcoming experiments with de- 
generate mixtures of 40 K and 87 Rb. This work has been sup- 
ported by the ESF under project BEC2000+, the INFM, the 
INFN, and the Italian Ministery for Scientific Research, under 
project PRIN-COFIN 2002. We dedicate this work to Noam 
Chomsky on his 75th birthday, for his tireless struggle in fa- 
vor of free research, human dignity and justice throughout the 
world. 



[1] M. Greiner et al, Phys. Rev. Lett. 87, 160405 (2001). 

[2] C. Orzel era/., Science 291, 2386 (2001). 

[3] M. Greiner et al, Nature 415, 39 (2002); M. Greiner et at, 

Nature 419, 51 (2002). 
[4] G. Modugno et at, Phys. Rev. A 68, 01 1601(R) (2003). 
[5] D. Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998). 
[6] W. Zwerger, J. Opt. B: Quantum Semiclass. Opt. 5, S9 (2003). 
[7] W. Hofstetter et al, Phys. Rev. Lett. 89, 220407 (2002). 
[8] A. Albus, F. Illuminati, and J. Eisert, Phys. Rev. A 68, 023606 

(2003). 

[9] H. P. Buchler and G. Blatter, Phys. Rev. Lett. 91, 130404 
(2003). 

[10] M. Lewenstein et al, Phys. Rev. Lett. 92, 050401 (2004); H. 

Fehrmann et al, Opt. Express 12, 55 (2004). 
[11] R. Roth and K. Burnett, Phys. Rev. A 69, 021601(R) (2004). 



[12] M. Cramer, J. Eisert, and F. Illuminati, cond-mat/03 10705 

[13] M. P. A. Fisher et al, Phys. Rev. B 40, 546 (1989). 

[14] L. P. Pitaevskii and S. Stringari, Science 298, 2144 (2002). 

[15] K. M. O'Hara et at, Science 298, 2179 (2002). 

[16] L. Viverit, C. J. Pethick, and H. Smith, Phys. Rev. A 61, 053605 

(2000); L. Viverit, ibidem 66, 023605 (2002). 
[17] R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. 

Phys. 62, 113 (1990). 
[18] R. T. D. Scalettar et al, Phys. Rev. B 39, 4711 (1989). 
[19] F. Schreck etal, Phys. Rev. Lett. 87, 080403 (2001). 
[20] T. Loftus etal., Phys. Rev. Lett. 88, 173201 (2002). 
[21] G. Roati et al, Phys. Rev. Lett. 89, 150403 (2002). 
[22] A. Simoni et al, Phys. Rev. Lett. 90, 163202 (2003). 



